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In reference (1) San Andres employs the NBS software package MIPROPS to 
account for density's dependence on pressure in the simulation of liquid annular 
seals. His example on a LH2 seal showed a significant change in the mass c °efficie 
compared to a constant density model. San Andres Yang and Childs (2,3) extend 
this analysis by including the pressure and temperature dependence of density, 
specific heat, viscosity, volumetric expansion and thermal conductivity in a coup 
solution of the energy, momentum and continuity equations. Their exampk s owe 
very significant changes in stiffness and inertia for a high speed (38,000 rpm), large 
L/D ratio (0.5) LOX seal, as compared to their constant temperature results. 

The current research rederived the San Andres-Yang-Childs (SYC) analysis 
and extended it to include not only the Moody friction model of SYCbut also the 
Hir’s friction model. The derivation begins with obtaining the local differential 
equations of continuity, momentum and energy conservation in the seal. Th ese 
equations are averaged across the film thickness to obtam the resulting 'bulk flow 
differential equations. Shear stress and convective heat loss throughthe stator (seal) 
and rotor are related to the Moody and Hir’s friction factor model. The Holman 
analogy is employed to relate heat conduction in or out of the fluid film s boundary 
layer to the friction induced shear stress. 

The steady state problem (d/dt=0) was solved using a shooting algorithm for 
the two-point boundary value problem. This require a simultaneous integration of 
the two momentum equations and the continuity and energy equation. Theresulte 
for temperature increase through the seal shows excellent agreement with the SYC 
model results as shown in figure 1. The SYC papers also describes an approximate 
solution algorithm which assumes constant properties and friction factors along the 
length of a concentric, straight seal. This model was deciphered and programmed 
and shows excellent agreement with the published SYC approximate solution 
results, a comparison of which is shown in figure 1. 

The linearization coefficient expressions were derived to solve the first order 
(perturbation) problem for the dynamic coefficients. This linearization procedure 
was performed for both the Hir’s and Moody models and revealed two errors in the 
SYC linearization coefficients for the viscosity and density in the circumferential 
momentum equations, and a missing convective heat flux term in the energy 
equation. The results showed that the Hir’s model linearization coefficients were 
quite different from their Moody counterparts, while maintaining a similar form as 
regards to programming. 

The non-dimensional equations employed in the preceding analysis were 
used to derive similarity conditions and expressions to infer LOX seal characteristics 
from those of a similar water seal. The branch is currently developing tius tester and 
required sizing information along with equations which reiate characteristics of the 
two seals. The similarity analysis was confirmed by running the TAMUSEAL code 
for a LOX seal and for its "similar" water seal. The results of these two runs showed 
nearly perfect agreement with those predicted by the similarity equations. This 
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numerical check was performed for both a Hit's and a Moody model type seal. The 
same study identified non-dimensional dynamic coefficients which remain invariant 
for seals that are mutually similar, i.e., obey the same conditions of similarity. 

The detailed analysis and results of this work may be found in the 430 page 
report, "Thermal and Similarity Studies for Cryogenic Liquid Annular Seals" issued 
by the Summer Faculty Fellow to the Mechanical Systems Controls Branch»Future 
work includes programming the first order solution to the thermohydrodynamic 
problem to obtain the resulting dynamic coefficients, including seal housing 
flexibility and extending the bulk flow model to include impeller forces. 

The Fellow also planned an installation of an impact damper on the TTB- 
ATD-HPOTP. The proposed location of the impact damper is shown in figure 2. 
This device will consist of 12-20 specially designed, cylindrical impactors contained 
in a ring type fixture. This type of damper has been successfully employed in LN2 
at Texas A&M.Testing of the impact damper may begin as early as Summer '94 if 
approved by the TTB Review Panel. 
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Figure 1 - Comparison Between the Exact and Approximate 
Temperature Rises 



Figure 2 - Proposed Location of the SSME-ATD-HPOTP 
Impact Damper 
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